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Abstract: The thermodynamics expected of systems undergoing third order phase transition 
has been investigated by identifying the orders through the analytic continuation of the 
functional of the free energy, using Ehrenfest thermodynamic theory. We developed the 
Euler - Lagrange equations for the order parameter and the vector potential and solved them 
for the first time using well - known mathematical formulations. 

Introduction 

The general theory of the thermodynamics in materials, especially systems undergoing phase 
transition(s) are generally complicated. Mean field theories have provided a way to analyze 
these systems, and has remained a vital tool in studying complex systems like the 
superconductors [1,2]. Even at these opportunities provided by mean field theories, for 
systems undergoing higher order phase transitions (>2), there abound many divergent views 
in literature to the existence and non-existence of phase transition of orders strictly greater 
than II [see 3,4,5]. Certainly, third order phase transition exists (even higher orders) since, for 
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a phase transition to be a second order, the relation j[SC e (H,T}/T~^iT = must hold. 

o 

However, as observed in many superconductors (for instance Na x CoOr yH 2 and BaKBiOi) 
[6,7,8,9,10,11,12]; Lattice and continuum gauge theories [13]; Chiral model [14]; Spin 
glasses [15]; DNA under mechanical strain [16]; Ferromagnetic and antiferromagnetic spin 
models with temperature driven transitions [17]; Invar types alloys [18]; and many other 
materials and systems, this conservation law is violated. Thus, we assert that, the non- 
detection of phase transitions of orders strictly greater than two might have been due to the 
erroneously hasty generalization that all that departs from phase transition of order two can 
always be explained in terms of thermodynamic fluctuation. 

The aim of present study is thus, to develop the thermodynamics expected of systems 
undergoing third order phase transition and most importantly, to solve the highly non - linear 
Euler - Lagrange equations for both the order parameter and the vector field using the system 
parameter from the local free energy of the functional of the free energy. The solution of 
these Euler - Lagrange equations has not been solved before by any group to the best of our 
knowledge. 

Theoretical Framework 

We will adopt a theoretical approach based on the classification of order of phase 
transitions by Ehrenfest [1 1,18,19] which is in terms of the ordered state free energy near the 
phase boundary and the functional of the system in terms of the local free energy. The 
validity of using Ehrenfest classification of order of phase transition has been questioned by 
many authors [20,21] which has erroneously led to the classification of order of phase 
transitions into first order and continuous phase transitions. However, as has been shown by 
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Hilfer [22], that rewriting the singular part of the local free energy within a restricted curve 
through the critical point in terms of the finite difference quotient and analytic continuation 
in p allows one to classify continuous phase transitions precisely according to their orders. 

According to classification of order of phase transitions by Ehrenfest, the transition in 
general can be of any order [11,22]. In a second order phase transition, the specific heat and 
the susceptibility, which are the second derivatives of the free energy (with respect to 
temperature and magnetic field, respectively) are discontinuous at the transition line [23]. 
Often, this discontinuity is replaced by a weak (or in general logarithmic singularity). Thus, 
one could view the Ehrenfest definition of an order of phase transition as one where the 
lower derivatives of the free energy are continuous at the transition but the higher derivatives 
are discontinuous. For a typical third order phase transition, all second order derivatives are 
continuous and all fourth order derivatives are singular at the transition. The third order 
derivatives are either discontinuous or a weakly logarithmic singular. 

The equilibrium thermodynamics is completely determined by the function, F\_T,(y/)~\ 

where (y/) is the local order parameter [24,25]. However, F must be invariant under the 

symmetry group say, G of the disordered phase in order to minimize the total energy [1]. 
Examples of such behavior are commonly found in magnetic domains in ferromagnetic [18] 
where the Ising (N = 1) Hamiltonian has been applied [21]. In general, F is a very 
complicated functional of (y/) but since in the vicinity of critical line, (y/) is essentially zero 

for T > T c , we will follow William - Bragg' s theory [24] and expand F in a Taylor series in 

y/). To make (y) to be spatially uniform in equilibrium in the ordered phase, we 

essentially for all cases redefine the order parameter. This suggests that F be expressed in 
terms of a local free energy density, /[r,(^(x))] which is a function of the field (w{x)} at 
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the point x only and also a part producing the energy cost for deviations from spatial 
uniformity [24]. After coarse graining [1,24], in its simplest form based on lattice 
formulation of Ginzburg - Landau theory [1,26], F is give by 

F = \d d xf[T,{yy{x))] + l -\d d xc{y{ W {x))) 2 (1) 

where c is a phenomenological constants known as the coupling constant or the gradient term 
and /[r,^(x))] is the local free energy density, which is well known in the y/ 4 - theory . It 

must be asserted that Eq. (1) fails in the short X - limit. However, we make haste to add that 
the short X - limit is never a problem as most system undergoing higher order phase 
transitions in the vicinity of the critical points are controlled by long - X fluctuations and as 
such, Eq.(l) is adequate for our problem so long as the short - X is significantly suppressed 
[27,28]. This is accomplished by the so - called hard cutoff mechanism where excitation of 
wave - number greater than a cutoff A □ , (where a is a length of order of the same 
range as the inter-particle interactions) are restricted in the partition trace [1]. 
Eq. (1) will be generalized to suit the thermodynamic theories of higher (>2) order phase 
transitions that will be developed, which is totally different from those developed for second 
order. 

Following [19,29,30] and using Eq. (1), the free energy can be generalized to Landau - 
like from as: 

F p (^r)=|jVH 2( ^ 2) {-vH 2+ ^H 4+ sl v H 2 } v ^> 2 ( 2 ) 

where the integer p is the order of the phase transition, a p = a" p (l-T/T c ) , bU , the term 
(VxA) 2 has been removed for brevity and it should be noted that due to gauge 
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transformation V-tV + lmA/fa ; where <f> =h/2e is the superconducting flux quantum and^ 
is the magnitude of the vector potential. Note that F is a functional and as such, the 
corresponding (Gibb's) free energy is to be identified with its minimum. For p = 3, Eq.(2) 
becomes: 



F 3 (y/ = Ae Up ) = \d d r 



II 4 7 I I 6 

-a 3 tm\ + b 3 \y/\ +c 3 



1 



2»o 



(Vxi) 2 



(3) 



Eq.(3) is the functional of the local free energy that contains all the thermodynamical 
information of the system. 

Analysis and Solution to the Functional of the System 

A critical look at Eq.(3) shows that it contains two important parameters that 
describes the thermodynamics of the system: the order parameter, y/ and the vector potential, 
A. We will solve the local free energy of Eq.(3) using variational principles: 1) with respect to 
the order parameter and 2) with respect to the vector potential respectively with and without 
any boundary conditions. 

Case I: Neglecting Gradient Term and Absence of any Field 

In the absence of any field and neglecting the gradient (coupling) term, we can write the local 
free energy of Eq.(3) as: 



f 3 =-a 3 t\y/\ 4 +b 3 \y/\ 6 
which, on minimization reduces to 



(4) 
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Using Eq.(5) in (4), 
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In its condensed degree of freedom, we can investigate the various thermodynamic quantities 
of interest using Eq.(6). The entropy, S-d^F^/dT and the specific heat, 

C — —T3 2 (F) jdT 2 are respectively 



_d(F)__4 al 



dT 



9 b 2 T c 



and 
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The order parameter, y/ can be expressed in terms of the local free energy (see Eq.(6)) as 



f( T ) = jV 6 o ; V?,=2a 3 t/3b 3 ; t >0 (t = 0; t <0) 



(8) 



Case II: Presence of Gradient Term and Field 

Next, we investigate, the effect of presence of filed and the incorporation of the 
gradient term to the thermodynamic behavior of the system. Again, with recourse to Eq.(3), 
we apply variational principles first with respect to the order parameter and secondly with 
respect to the potential field. 

In terms of the order parameter, we obtain that 



2 \^¥\ 3 

c 3 V y/ + c 3 - - + 2a 3 ty/ -3b 3 y/ = 

¥ 



(9) 
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which is the Euler - Lagrange equation of the free energy for the order parameter: y/ = Ae" p 
in the third order phase transition regime. Also, for the vector potential, A, we have, 



1 6* 
7T A = - 
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where e* and m are the effective charge and effective mass respectively. In the Meissner 
state at the depth X, the gradient term can be neglected to obtain the Euler - Lagrange 
equation for the vector potential in terms of the penetration length, A, as: 



r 2 = e 



*2 
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c m 



2n^ 



W\ 
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which, is the Euler - Lagrange equation in terms of the penetration length; and can be 
rewritten using gauge transformation as, 



V 2 A-2 Mo c 



2n^ 



\y/\ A = 



(12) 



which is nothing but the London equation for the behavior of the system in third order phase 
transition regime. 

The specific heat can be rewritten in terms of the penetration length by using Eq.(l 1) 
in Eq.(8) as 

f a V p3 



C 3 (T) = —W- 
1 ) 2(2 Mo c) 
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Solution to the Euler - Lagrange Equations 

We will seek the solutions to Eqs.(9 &12) which are highly non - linear Euler 
Lagrange equations for the order parameter and the vector potential respectively. 



As a starting point, Eq.(12) can be redefine in terms of the potential, ^(r) as 



V 2 A + n(r)A = 



(14) 



where Q = -2 ju c 
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Recall that V 2 = — - + — - + — - , and the function Q(V) will be regarded in this case as a 
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potential function independent of and <j> when we transform the Cartesian coordinates 
(x,y,z) to polar (r,0,(f) . Thus, Eq.(14) reduces to 
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which on using the method of separation of variables: A - R(r}S(0,(f) , gives 
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Since the LHS is depends only on r and RHS is independent of r , therefore each side must 
be equal to the same constant say, // for the equation to be valid. Thus, from Eq.(17), we 
have for the LHS that; 
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Provided Q(r) - X 2 , Eq.(18) reduces to Bessel's equation, with general solutions; 



R (r) = A,J M (/T 2 r ) + B X J_ M (/T 2 r ) for fractional ju 



(19a) 



R(r) = A X J M (/T 2 r) + BJ^ (/T 2 r ) for integral ju 



(19b) 



Similarly, taking the RHS of Eq.(17), 
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Again, applying separable, S(0,(f>) - ©(^)O(^), Eq.(20) reduces to 
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where / is constant of separation. Taking the RHS of Eq.(21), we have; 
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with the solution 



0(^) = e i¥ =cos/^+/sin/^. 



(23) 



For the LHS of Eq.(21), we have; 
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For convenience we have set /u 2 = n(n + 1) . Making change of variable; i.e. cos# = u so that 
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, then Eq.(24) reduces to; 
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Eq.(25) is the associated Legendre equation, and has its solution as; 
= AP; {u)+BQ: (u) = AP; (cos0)+ BQ: (cos0) 



(26) 



where A and B are constants to be determined by the initial conditions. Hence, we obtain 
the complete general solution to the Euler - Lagrange equation for the vector potential as 

A{r,6,(j)) = (A { J M (Xr) + B x J_ M (/Lr))(cos/^ + /sin/^)(j 2 P B '" (cos 0) + B 2 Q™ (cos&)) (27) 

For the order parameter, since c is a gradient term, we can rewrite Eq.(9) as: 
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Solving in ID, 
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dx y/ 
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= 0, y/ = y/(t,x) 



(29) 



Assuming a constant temperature gradient field (such that Hsa constant) where the 
variation of the system depends entirely on the order parameter, further assuming natural 
units, c = 1 and also noting that b >1 and as such also a constant, we can rewrite Eq.(29) as: 
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d 2 y/ if dy/ 
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V dx j 



+ 2y/-3y/ 3 =0 (30) 



To linearizing the non - linear differential (Eq.(30)) we use a numerical procedure by letting 
Zj = y/ , thus reducing to, 



z"+ — [z'J +2 Zl -3 Zl 3 =0 (31) 



Next, let let z 2 = z[ hence we have 
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Observe that Eq.(28) has been reduced to a system of first order differential equations which 
can then be solved. 

Conclusion 

The thermodynamics of third order phase transition have been analyzed using mean 
field theory based on Ehrenfest thermodynamics. Certain thermodynamical relations well - 
known for systems undergoing second order phase transitions are recovered for third order 
phase transition. We developed equations for the phase boundary in terms of discontinuities 
in thermodynamic observables, derived Ginzburg - Landau free energies, developed Euler - 
Lagrange equations for the order parameter and the vector potential, and solved the highly 
non - linear Euler - Lagrange equations for both the order parameter and the vector potential. 
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